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(57) In a wheel condition estimating apparatus for 

estimating a wheel condition in a wheel resonance sys- 
tem including a frictional characteristic between a tire 

and a road surface, a transfer characteristic of the wheel 

resonance system from the vibration input to the output 

response is expressed by a transfer function including, 

as an unknown component of a wheel condition, a phys- 
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ical quantity relating to ease of a slipping between the 
tire and the road surface, an output response with 
respect to the vibration input to the wheel resonance 
system is detected and the unknown component which 
substantially satisfies the detected output response is 
estimated on the basis of the transfer function. 
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Description 

BACKGROUND OF THE INVENTION 
5 Field of the Invention 

The present invention relates to a wheel condition estimating apparatus, and more particularly to a wheel condition 
estimating apparatus for estimating a wheel condition from an output response or the like of a wheel resonance system 
on the basis of a vibration model of the wheel resonance system including a f rictional characteristic between a tire and 
10 a road surface. 

Description of the Related Art 

In recent years, in accordance with an increase in efforts to improve the safety of automotive vehicles, a preventive 
15 safety techniques have been researched and developed. An antilock brake systems (ABS), which are a representative 
safety apparatus, have already been used in many vehicles. 

ABS control on the basis of a new principle which pays attention to the resonance phenomenon of the tire has been 
proposed, and the study thereof continues (Japanese Patent Application No. 7-220920 and the like). This technique is 
such that a micro vibration having a frequency component equal to a resonance frequency when the tire grips a road is 
20 applied to the brake pressure, and the average brake pressure is controlled on the basis of a resonance gain (micro 
amplitude of wheel speed at the resonance frequency/vibration amplitude of brake pressure) of a tire resonance system 
at that time. 

It is known that the resonance gain is a physical quantity related to the gradient of a coefficient of friction \i with 
respect to a slip ratio S (hereinafter referred to as a n gradient) of a coefficient \i of friction in a so-called S-n character- 
's istic (a changing curve of a coefficient of friction n with respect to a Slip ratio S). It is expected that a frictional charac- 
teristic, which relates to ease of slipping between a tire and a road surface at a time of braking, can be estimated on the 
basis of the resonance gain. 

However, since the prior art mentioned above presupposes that the resonance frequency of the tire resonance sys- 
tem is known, there is a problem that it cannot immediately respond in the case that the resonance frequency changes 
30 due to the tire being replaced or the air pressure in the tire changing. For example, when the air pressure in the tire falls, 
the resonance frequency at the time when the tire grips the road surface changes to a lower side. Therefore, the value 
of the resonance gain becomes small, and the accuracy of detecting whether the current state is a state immediately 
before a peak n value by comparison with a predetermined reference gain deteriorates. 

Further, in the above-described prior art, the brake pressure is minutely vibrated at a time of calculating the reso- 
35 nance gain. However, because the micro vibration can be applied only at a time of braking due to the structure of a 
brake portion, there is a problem that the opportunity to estimate the road surface condition is limited to the time of brak- 
ing in the case of a vehicle having the brake portion of this type. Still further, there are cases in which micro vibration of 
the brake pressure is not preferable such as a time of driving, a time of steady running and the like. 

40 SUMMARY OF THE INVENTION 

The present invention has been achieved by taking the facts mentioned above into consideration, and an object of 
the present invention is to provide a wheel condition estimating apparatus which can estimate a road surface condition 
equally in various running conditions such as braking, driving, steady running and the like, and which can always esti- 

45 mate with high precision a frictional characteristic which relates to ease of slipping between a tire and a road surface 
regardless of variations in a resonance frequency due to replacement of a tire or changes in air pressure of a tire, and 
which enables highly-accurate ABS control on the basis of the frictional characteristic. 

In order to realize the aforementioned object, in accordance with a first aspect of a present invention, there is pro- 
vided a wheel condition estimating apparatus for estimating a wheel condition in a wheel resonance system including 

so a frictional characteristic between a tire and a road surface, comprising: detecting means for detecting an output 
response with respect to a vibration input to the wheel resonance system; and estimating means for expressing a trans- 
fer characteristic of the wheel resonance system from the vibration input to the output response by a vibration model 
including, at least a physical quantity relating to ease of a slipping between the tire and the road surface as an unknown 
component of a wheel condition, and for estimating at least the unknown component which substantially satisfies the 

55 output response detected by the detecting means on the basis of the vibration model. 

Further, in accordance with a second aspect of the present invention, the wheel condition estimating apparatus of 
the first aspect further comprises modifying means for modifying the vibration model on the basis of the output 
response detected by the detecting means. 
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Still further, in accordance with a third aspect of the present invention, in the wheel condition estimating apparatus 
of the first aspect or the second aspect, the vibration model further includes, a physical quantity relating to a constant 
of a torsion spring of the tire as an unknown component of the wheel condition. 

Furthermore, in accordance with a fourth aspect of the present invention, the wheel condition estimating apparatus 
s of any of the first aspect to the third aspect further comprises vibrating means for applying a vibration input to the wheel 
resonance system. 

Moreover, in accordance with a fifth aspect of the present invention, the wheel condition estimating apparatus of 
the fourth aspect further comprises vibration input detecting means for detecting the vibration input applied to the wheel 
resonance system by the vibrating means. 

10 Further, in accordance with a sixth aspect of the present invention, the wheel condition estimating apparatus of any 
of the first to the fifth aspects further comprises calculating means for calculating a determining value for determining 
whether or not the output response detected by the detecting means is a periodical output on the basis of the output 
response, and selecting means for selecting the output to be delivered to the estimating means among the output 
responses detected by the detecting means on the basis of the determining value calculated by the calculating means. 

15 wherein the estimating means uses the output selected by the selecting means at a time of estimating the unknown 
component. 

Still further, in accordance with a seventh aspect of the present invention, in the wheel condition estimating appa- 
ratus of the sixth aspect, the calculating means calculates a correlation coefficient between the output response 
detected by the detecting means and a converted value obtained by converting the output response such that the peri- 
20 odical portion of the output response changes in correspondence to the change of the output response, as the deter- 
mining value. 

(Principles of the Present Invention) 

25 Principles of the present invention will be described below with reference to Figs. 1 to 3. Here, Fig. 1 shows an 
equivalent dynamic model of a wheel resonance system, Fig. 2 shows a frictional characteristic between a tire and a 
road surface, which defines a transfer characteristic of the wheel resonance system shown in Fig. 1 , and Fig. 3 shows 
an example of a vibration model from a vibration input to an output response in the transfer characteristic of the wheel 
resonance system shown in Fig. 1 . 1 

30 First, as shown in Fig. 1 . a vibration phenomenon in a wheel when a vehicle runs at a vehicle speed V (an angular 
velocity gv), that is. a vibration phenomenon of a wheel resonance system formed by at least a wheel and a road sur- 
face, is considered by referring to a dynamic model which is obtained by equivalent modeling by using a wheel rota- 
tional axis. The various quantities shown in Fig. 1 are as follows. 

35 J-j : moment of inertia at rim 
J 2 : moment of inertia at belt 
K: constant of torsion spring of tire 

T 1 : braking and driving torque (driving torque is positive, braking torque is negative) 
<o 1 : angular velocity at rim 

40 cog: angular velocity at belt 

e s : torsion angle between rim and belt 

T d : disturbance from road surface 

T L : force generated between tire and road surface 

In the dynamic model of a wheel resonance system in Fig. 1 , the brake and drive torque T-, applied to the rim 

45 is transmitted to the belt through the torsion spring constant K of the tire, and is applied to the road surface through 
the belt surface. At this time, the generated force T L serving as the reaction of the brake and drive torque is 
applied to the wheel from the road surface, with a contact point between the belt and the road surface being a ref- 
erence point. 

so This generated force T u is caused by the frictional force between the tire and the road surface, and acts in a direc- 
tion opposite to the direction of the brake and drive torque T v Accordingly, the generated force T L acts in a direction 
opposite to the rotational direction of the wheel (the direction of c^) in the case that the drive torque T t is applied to the 
rim at a time of driving, and the generated force T L acts in the rotational direction of the wheel in the case that the brak- 
ing torque T t is applied at a time of braking. Further, if the road surface is uneven, a torque of a road surface distur- 

55 bance AT d generated by the unevenness also acts on the tire. 

In this case, when the brake is applied at a time when the vehicle runs at a certain speed V (the value when con- 
verted into a rotation system is gv), a slip is generated between the tire and the road surface. However, the generated 
force T L generated between the tire and the road surface at this time changes like the functional relation shown in Fig. 
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2 with respect to the slip ratio S 1 expressed by the following formula (the slip ratio is in a positive area). 

^ = " (1-1) 



10 Similarly, also in the case that the driver steps on the acceleration pedal so as to accelerate at a time when the vehi- 
cle runs at a certain speed V, a slip is generated between the tire and the road surface. However, the generated force 
T L at this time changes like the functional relation shown in Fig. 2 with respect to a Slip ratio expressed by the fol- 
lowing formula. 




(1-2) 



20 

In this case, when the rotational direction of the wheel is the forward direction, the generated force T L between the 
tire and the road surface can be expressed by the following formulae. 

25 At the time of braking: T L = WFV(S 1 ) (2) 

At the time of driving: TL = -WRja(S 2 ) (3) 

Here, W is the wheel load, R is the radius of the dynamic load of the tire, and u is the coefficient of friction between 
30 the tire and the road surface. In this case, \i is expressed as a function of the slip ratio S 1 or S 2 . 

As shown by the S-n curve in Fig. 2, it can be understood that a relation is established in which the generated force 
T L is 0 when the slip ratio is 0, the generated force T L at the time of braking has a positive peak value at a certain pos- 
itive Slip ratio, and the generated force T L at the time of driving has a negative peak value at a certain negative slip ratio. 
Further, at various points of operation, the gradient of the generated force with respect to the slip ratio becomes a char- 
ts acteristic value, fa example, a value near 0 at a time of the peak value. Therefore, the ease of slipping between the tire 
and the road surface can be expressed by using the gradient. 

When the brake and drive torque applied to the rim is vibrated around the average brake and drive torque T 1 by an 
amplitude AT 1 in the dynamic model in Fig. 1, the vibration torque component is present as the vibration component 
Ag>i around the wheel speed Further, in the case that the vibration component AT d is present in the road surface 
40 disturbance T d . the vibration component generated by the disturbance is also added to the vibration component Ao-j of 
the wheel speed. 

When the transfer characteristic of the wheel resonance system in Fig. 1 is expressed by the vibration model in the 
various points of operation in Fig. 2, the following formula can be established. 

45 Aio 1 = H 1 (s) AT 1 + H 2 (s) AT d (4) 

Here, the following formulae can be established. 



50 



55 



H 2 (s) = 



J?s 2 +D„S+K 



K 

J X J^ + J x DoS 2 + K (J x + J 2 )s + KD 0 



(5) 
(6) 
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Here, s denotes a Laplace operator. 

Further, D 0 is expressed by D 10 and D 2 o in the following formulae, in accordance with the time of braking and the 
time of driving, respectively. 



D 10 - WR 



CO 



10 



ds, 



2 <?j" 



(7) 
(8) 



15 



S 10 and S 2 o represent Slip ratios at a certain point of operation at a time of braking and at a time of driving, respec- 
tively, and cvo denotes a vehicle speed at that point of operation. 

In accordance With formula (7), D 10 is proportional to the wheel load W and to the Gradient (d\\J dS^) of the S-u. 
20 curve at the Slip ratio S 10 at the point of operation, and is inversely proportional to the vehicle speed oovo at the point of 
operation. Further, when S 20 >s near 0, the same also holds for D 20 . 

It is supposed that depends upon the Slip ratio in the above description. However, in the case that \i depends 
upon the slip speed, the following formulae can be established by redefining St = cd v - g> 2 and S 2 = a> 2 - <o v . 



25 



30 



D l0 - WR 



D 2Q - WR 



(9) 



(10) 



35 In this case, D 10 and D 20 are proportional to the wheel load W and to the Gradient of the S-ji curve at the Slip ratios 
S 10 and S 2 o at the operating point. 

Since the above-described vibration model expresses the operation at an arbitrary point of operation, it also 
describes a special case of steady running in which no braking and nor driving is performed. In the steady running, the 
point of operation becomes an origin of the S-^i curve, and D 0 = D 10 = D 20 expresses the n Gradient at the origin. 

40 Further, the above-described vibration model is expressed by a third-degree system with respect to a Laplace oper- 
ator s. However, in order to describe the physical phenomenon of vibration, a second-degree system is sufficient. When 
the third-degree model is approximated to the second-degree model, the following formulae can be obtained. 



45 
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H 2 (s) = 



JfitS 2 + K (jr, + J 2 )s + KD Q 

K 

JiDjS 2 + K (J, + J 2 )s + KD 0 



(11) 
(12) 



55 It can thus be understood that the above-described vibration model expresses a response of the wheel speed 
vibration (Ao^) as the output response with respect to the vibration input torque (AT.,) to the resonance system and the 
road surface vibration (AT d ) generated to the wheel resonance system by the tire rotating on an uneven road surface, 
in the wheel resonance system including the frictional characteristic between the tire and the road surface. It can also 



5 



<EP 0891904A2J_> 



< J 

EP 0 891 904 A2 

be understood that the vibration model includes a physical quantity D 0 relating to the ease of slip between the tire and 
the road surface. 

Fig. 3 illustrates the state of the transfer from the vibration input to the output response in the above-described 
vibration model. 

5 The appropriateness of this vibration model can be shown by experimental results. Figs. 4A and 4B show results 

obtained by conducting an experiment on the frequency characteristic from the vibration input to the wheel speed vibra- 
tion at the time when, at a given constant vehicle speed, the vibration input is added to the brake pressure P m so as to 
effect braking, with respect to various values (0.98 [MPa] to 4.90 [MPa]) of the brake pressure P m . Fig. 4A shows the 
frequency characteristic (amplitude characteristic) of the ratio of the wheel speed vibration with respect to the vibration 

10 amplitude of the brake pressure P m (resonance gain), and Fig. 4B shows the phase characteristic from the vibration of 
the brake pressure P m to the wheel speed vibration. 

It can be understood from Fig. 4A that the peak value of the resonance gain near the resonance frequency (about 
40 Hz) decreases as the brake pressure increases. Further, it can be understood from Fig. 4B that in the phase from 
the brake pressure to the wheel speed vibration, the phase characteristic in the case of setting a value near the reso- 

15 nance frequency to be a boundary differs greatly in accordance with the magnitude of the brake pressure P m . 

The slip ratio and the braking force are increased and thus the gradient of the S-\i curve is reduced by an increase 
in the brake pressure P m . Therefore, the increase in P m and the decrease in the gradient D 0 of the S-ji curve have sub- 
stantially the same physical meaning. Accordingly, even if Dq is varied, it can be expected that results similar to those 
of Figs. 4A and 4B can be obtained. 

20 On the other hand, when the response characteristic of the wheel speed vibration with respect to the brake vibra- 
tion input is calculated on the basis of the above-described vibration model of the present invention, the amplitude char- 
acteristic shown in Fig. 5A and the phase characteristic shown in Fig. 5B are obtained. It can be understood from Fig. 
5A that when the Gradient D 0 of the S-n curve is reduced at the resonance frequency (40 Hz), the resonance peak is 
reduced in the same manner as shown in Fig. 4A. Further, it can be understood from Fig. 5B that the phase character- 

25 istic in the case of setting a value near the resonance frequency to be a boundary differs greatly in accordance with the 
magnitude of the Gradient D 0 . so that there is a characteristic similar to Fig. 4B. Accordingly, the vibration model of the 
present invention well expresses the transfer characteristic of an actual wheel resonance system. 

Further, Figs. 6A and 6B show frequency spectra of the wheel speed vibration of drive wheels (a right rear wheel 
and a left rear wheel) obtained at a time of actual running on an asphalt road (a dry road) and a dirt road, respectively. 

30 It can be understood that in the case of an asphalt road (Fig. 6A), the resonance peak clearly appears. However, 
no resonance peak appears in the case of a dirt road (Fig. 6B) on which the tire runs in a state near racing (idle rotation) 
due to power slide and the like. 

Figs. 7A and 7B respectively show the amplitude characteristic and the phase characteristic of the wheel speed 
vibration with respect to the tire road surface vibration input on the basis of the vibration model of the present invention. 

35 It can also be understood from these drawings that since the resonance peak comes to disappear as the gradient D 0 
of the S-\i curve decreases, the amplitude characteristic and the phase characteristic well coincide with the character- 
istics of the experimental results. 

As described above, the present invention is structured so as to estimate the physical quantity D 0 relating to the 
ease of slipping between the tire and the road surface, on the basis of the vibration model which well expresses the 

40 actual transfer characteristic. Accordingly, on the basis of the above-described vibration model which includes, as the 
unknown component of the wheel condition, at least a physical amount relating to the ease of slipping between the tire 
and the road surface, the estimating means of the present invention estimates at least an unknown component which 
can substantially satisfy the output response detected by the detecting means. 

Further, as in the third aspect, in the case that the vibration model further includes, as an unknown component of 

45 the wheel condition, a physical quantity relating to the torsion spring constant K of the tire, the torsion spring constant 
K of the tire and the physical quantity relating to the constant can be estimated by the estimating means. The resonance 
frequency of the wheel resonance system and the like may be used as the physical quantity relating to the torsion spring 
constant of the tire. 

Next, the principles of estimation by the estimating means in accordance with the present invention will be 
so described. 

As shown in Fig. 3, in the case of expressing the vibration model of the present invention by a transfer function 1 
and a transfer function 2, the estimation of the unknown component by the estimating means is equivalent to the iden- 
tification of these transfer functions. 

A z-transformation discrete model is often employed as a transfer function for identification. However, there are the 
55 following problems in the case of identifying a discrete time model. 

(1) Since the identification accuracy of the discrete time model is dependent on the sampling period, many trials 
and errors are required for obtaining a suitable sampling period. 
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(2) After identifying by the discrete time model, it is necessary to inversely convert to a continuous time model in 
order to calculate the physical quantity forming the system. However, a higher function is required for the calcula- 
tion, so that the calculating time and the calculating errors are increased. Further, the inverse conversion is not uni- 
formly obtained. 

5 

Accordingly, the z-transformation discrete model is not suitable to a model for calculating the physical quantity cor- 
responding to the ease of slipping on the road surface. 

Therefore, in the present invention, a continuous time model is identified. In the identification of the continuous time 
model, there are none of the problems caused by the identification of the discrete time model, so that there is an advan- 
10 tage in that the physical quantity corresponding to the ease of slipping on the road surface can be directly calculated. 

For example, in the case of identifying the transfer function of a second-degree (quadratic) continuous time model, 
the following transfer functions obtained by modifying formulas (11) and (12) may be used as transfer functions to be 
identified. 

15 

(s) ^ b fl s'+b, s +b 2 (11 . 2) 

1 s a + + a. 



25 

At this time, as in the fourth aspect, in the case of giving the vibration torque AT 1 as the vibration input to the wheel 
resonance system, Aco 1 = AG 1 (s) AT 1 can be obtained by formula (4) when the road surface disturbance is ignored 
since the road surface disturbance is minute in comparison with the vibration torque. 

For example, in the case of using a least square method, the unknown parameter can be estimated by applying the 
30 least square method to each of the data obtained by successively applying the detected Acoi to a formula which modi- 
fies the above formula with respect to the unknown parameter [a^ a 2 ] or [a 1 a 2 b 0 b-i b 2 ] in accordance with a linear func- 
tion form. 

In this case, as in the fifth aspect, in the case in which the vibration input detecting means can detect the vibration 
input AT-, applied to the wheel resonance system by the vibrating means, all of the unknown parameters [a-j a 2 b 0 bj 
35 b 2 ] can be estimated. In contrast, in the case of not detecting AT 1f [a-i a 2 ] can be estimated. 

In this way, in the case of applying the vibration input to the wheel resonance system, the unknown parameter can 
be accurately estimated even when the tire is running on a good road surface which is for the most part even. 

In contrast, in a case in which the vibration torque A^ is not applied, the unknown parameter can be estimated by 
applying the least square method to each of the data obtained by successively applying the detected Ao 1 to a formula, 
40 which is obtained by modifying the formula Aoo 2 = AG 2 (s) AT d obtained by the formula (4) in the case of A^ = 0, to a 
linear function form with respect to the unknown parameter [a-\ a 2 ]. In this case, there is an advantage that the unknown 
parameter can be estimated even in the case that the vibration torque A^ cannot be applied. 

Any of various kinds of modified least square methods can be used for suppressing the identification error. The 
method of instrument variable, the enlarged least square method and the generalized least square method which have 
45 been conventionally known can be employed as the modified least square method. 

The estimating means in accordance with the present invention can estimate, as the following formula (13), the 
physical quantity relating to the road surface ix gradient D 0 from the corresponding relation between the transfer func- 
tions (1 1) and (11-2) or the functions (12) and (12-2) of the vibration model by using the estimated parameters a 1 and 
a 2 . 

50 



a, J x + J a 

55 

The estimating means can estimate the physical quantity relating to the torsion spring constant of the tire as a 2 = K/J , 
Further, the estimating means can thereby estimate the resonance frequency of the wheel resonance system. 
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In this way, when the physical quantity related to the \i Gradient D 0 of the road surface can be calculated, it can be 
judged that the frictional characteristic between the tire and the road surface is in a saturated state when the physical 
quantity is small, so that the ease of slipping on the road surface can be judged immediately. Further, the n Gradient D 0 
can be calculated without presupposing that the torsion spring constant of the tire is known. Therefore, the physical 

5 quantity relating to the ease of slipping between the tire and the road surface can be accurately determined regardless 
of whether the resonance frequency is changed due to replacement of the tire or a change in the tire air pressure. Fur- 
ther, the air pressure of the tire can be judged on the basis of the estimated value of the physical quantity relating to the 
torsion spring constant of the tire. 

Further, in accordance with the second aspect, the vibration model is modified on the basis of the detected output 

10 response. For example, in order to improve the identification accuracy, in the case that a pre-treating means having a 
frequency characteristic corresponding to the resonance frequency of the wheel resonance system is provided at a step 
before the transfer function of the vibration model, the parameter of the pre-treating means is properly changed together 
with calculation being carried out by the least square method on the basis of the detected output response, in this case, 
the parameter of the pre-treating means can be expressed by the parameter of the transfer function. The parameter of 

15 the transfer function can again be updated by the output response passing through the changed pre-treating means 
and the transfer function. In this way, since the frequency characteristic of the pre-treating means approaches the fre- 
quency characteristic which is appropriate to the resonance frequency corresponding to the change in the torsion 
spring constant of the tire, estimation can be performed with higher accuracy than that of a case in which the pre-treat- 
ing parameter is fixed. 

20 The calculating means in accordance with the sixth aspect of the present invention calculates the determining 
value for determining whether or not the output response is a periodical output on the basis of the output response 
detected by the detecting means. In this case, the calculating means may be structured such as to calculate the rela- 
tional coefficient between the output response detected by the detecting means and a converted value obtained by con- 
verting the output response such that the periodical portion of the output response changes in correspondence to the 

25 change of the output response, as the determining value. In this case, the conversion mentioned above is performed by 
any of an even number of differentiation, an even number of difference, an even number of integration and an even 
number of summation. 

In this case, even when the output response having an unexpectedly large magnitude is converted in the above 
manner, the converted value includes a lot of components except the periodical portions. Accordingly, since the value 

30 obtained by converting the unexpected output response in the manner mentioned above does not change in corre- 
spondence to the change of the output response, there is a small relationship between the unexpected output response 
and the value obtained by converting the output response in the manner mentioned above. 

Further, the value obtained by converting the output in which the periodical change does not significantly appear in 
the manner mentioned above has a small number of components converted such as to change in correspondence to 

35 the change of the output response. Accordingly, there is a small relationship between the output in which the periodical 
change does not significantly appear and the value obtained by converting the output response in the manner men- 
tioned above. 

In this case, the calculating means may be structured such as to calculate the determining value mentioned above 
at every time when a predetermined number of (for example, one) output response is detected by the detecting means. 
40 Then, the selecting means selects output responses to be inputted to the estimating means among the output 
responses detected by the detecting means on the basis of the determining value calculated by the calculating means. 

Accordingly, the estimating means uses the outputs selected by the selecting means at a time of estimating the 
unknown component mentioned above. 

As mentioned above, since the apparatus is structured such as to calculate the determining value for determining 
45 whether or not the output response is the periodical output, select output responses to be used for estimation among 
the detected output responses on the basis of the calculated determining value, and estimate the unknown component 
by using the selected outputs, the apparatus can estimate the unknown component mentioned above on the basis of 
the output response except the output response having an unexpectedly large magnitude and the output in which the 
periodical change does not significantly appear. 

50 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a schematic diagram which shows a dynamic model equivalent to a wheel resonance system in accord- 
ance with the present invention. 
55 Fig. 2 is a diagram which shows a frictional characteristic between a tire and a road surface in a wheel resonance 
system in accordance with the present invention. 

Fig. 3 is a schematic diagram of a vibration model in accordance with the present invention. 

Figs. 4A and 4B are diagrams which show results of an experiment on a frequency characteristic in a wheel reso- 
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nance system at a time of actually variously changing a brake pressure P m , wherein Fig. 4A shows an amplitude 
characteristic and Fig. 4B shows a phase characteristic. 

Figs. 5A and 5B are diagrams which show frequency characteristics of the vibration model in accordance with the 
present invention respectively calculated at times of variously changing a physical quantity D 0 relating to ease of 
5 slipping between a tire and a road surface in the vibration model, wherein Fig. 5A shows an amplitude characteristic 

and Fig. 5B shows a phase characteristic. 

Figs. 6A and 6B are diagrams which show results of analyzing a frequency of a wheel speed signal actually 
detected in a drive wheel, wherein Fig. 6A relates to an asphalt road and Fig. 6B relates to a dirt road. 
Figs. 7A and 7B are diagrams which show frequency characteristics of the vibration model in accordance with the 
w present invention respectively calculated at times of variously changing a physical quantity D 0 relating to ease of 
slip between a tire and a road surface in the vibration model, wherein Fig. 7A shows an amplitude characteristic 
and Fig. 7B shows a phase characteristic. 

Fig. 8 is a block diagram which shows a structure of a first aspect relating to a first embodiment. 
Fig. 9 is a block diagram which shows a structure of a second aspect relating to a first embodiment. 
15 Fig. 10 is a graph which shows results of estimating a physical quantity D 0 relating to ease of slipping between a 
tire and a road surface, which results are estimated by a wheel condition estimating apparatus in accordance with 
a first embodiment. 

Fig. 1 1 is a block diagram which shows a structure of a second embodiment in accordance with the present inven- 
tion. 

20 Fig. 12 is a graph which shows results of estimating a resonance frequency of a wheel resonance system, which 
results are estimated by a wheel condition estimating apparatus in accordance with the second embodiment. 
Fig. 13 is a block diagram which shows a structure of a first aspect relating to a third embodiment of the present 
invention. 

Figs. 14A to 14C are graphs which show a vibration input waveform of a vibration input means in accordance with 
25 the third embodiment, wherein Fig. 14A shows a pseudo random waveform, Fig. 14B shows a waveform in a step- 
ping manner, and Fig. 14C shows a waveform in an impulse manner. 

Fig. 15 is a block diagram which shows a second aspect relating to the third embodiment of the present invention. 
Fig. 1 6 is a block diagram which shows a first aspect relating to a fourth embodiment of the present invention. 
Fig. 1 7 is a block diagram which shows a second aspect relating to the fourth embodiment of the present invention. . 
30 Fig. 1 8 is a block diagram which shows a structure relating to a fifth embodiment of the present invention. 
Fig. 1 9 is a graph which shows two signals comprising a periodical signal and a non-periodical signal. 
Figs. 20A and 20C are graphs which show a non-periodical signal, and Fig. 20B is a graph which shows a period- 
ical signal. 

Fig. 21 is a graph which shows a relationship between an output response x and a value y obtained by converting 
35 the output response such that a periodical portion changes in correspondence to a change of the output response. 

Fig. 22 is a flow chart which shows a selecting operation routine performed by selecting means. 

Fig. 23A is a graph which shows a relational coefficient between a periodical signal and a non-periodical signal (a 

noise) in a time sequence manner, and Fig. 23B is a graph which shows a vehicle speed of the periodical signal 

and the non-periodical signal (the noise) in a time sequence manner. 
40 Figs. 24A and 24B are graphs which show a wheel speed signal of relatively low periodicity, and Figs. 24C and 24D 

are graphs which show a wheel speed signal of relatively high periodicity, wherein Figs. 24A and 24C relate to a 

front right wheel speed signal and Figs. 24B and 24D relate to a front left wheel speed signal. 

Fig. 25 is a graph which shows a relationship between frequency and spectrum of periodical signal and that of a 

signal in which periodical change does not appear remarkably. 
45 Fig. 26A is a graph showing, with passage of time, correlation coefficient of periodical signal and that of a signal in 

which periodical change does not appear remarkably. Fig. 26B is a graph showing, with passage of time, wheel 

speed corresponding to periodical signal and a signal in which periodical change does not appear remarkably. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

50 

Embodiments of a wheel condition estimating apparatus in accordance with the present invention will be described 
below in detail with reference to the drawings. 

(First Embodiment) 

55 

A wheel condition estimating apparatus in accordance with a first embodiment is structured such as to calculate a 
(j gradient in the case that only a road surface disturbance AT d is input to a wheel resonance system as a vibration 
input, and includes a first aspect and a second aspect. 
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A structure of a first aspect of the wheel condition estimating apparatus in accordance with the first embodiment is 
shown in Fig. 8. As shown in the same drawing, this wheel condition estimating apparatus is constituted by wheel speed 
defecting means 1 for detecting a wheel speed go-, of each of the wheels, a pre-f ilter 2 for detecting a wheel speed vibra- 
tion Aa> 1 of each of the wheels, as an output response of a wheel resonance system receiving a road surface distur- 
5 bance AT d from the detected wheel speed of each of the wheels, transfer function identifying means 3 for identifying 
the transfer function of each of the wheels which satisfies the detected wheel speed vibration Ao> 1 by using a least 
square method on the basis of the vibration model shown in Fig. 3, and u gradient calculating means 4 for calculating 
a gradient of a coefficient of friction between the tire and the road surface on the basis of the identified transfer func- 
tion. 

io in Fig. 8, the wheel speed detecting means 1 can be constituted by a so-called wheel speed sensor which outputs 
a sensor output signal in accordance with a wheel speed, and calculating means for calculating an actual rotational 
speed signal of each of the wheels from the sensor output signal. 

Further, the pre-f ilter 2 can be constituted by a band pass filter for passing only a frequency component at a fixed 
band area around a frequency expected as a resonance frequency of the present wheel resonance system, a high pass 
is filter for passing only a frequency component at a high band area including the resonance frequency component, or the 
like. In this case, in accordance with the first aspect, the parameter for defining the frequency characteristic of the band 
pass filter or the high pass filter is fixed to a constant value, and in a second aspect mentioned below, the parameter is 
changed in such a manner as to correspond to the parameter identified by the transfer function identifying means 3. 
In this case, an output of the pre-f ilter 2 is structured such as to remove a direct current component. Accordingly, 
20 only the wheel speed vibration Aoj around the wheel speed can be extracted. 
Here, a transfer function F(s) of the pre-filter 2 is set to be as follows. 



25 




(14) 



i-0 



30 Here, q is a coefficient of the filter transfer function, and s is a Laplace operator. 

Next, a calculating formula on which the transfer function identifying means 3 depends is introduced. In this case, 
in accordance with this embodiment, the operation is performed such that the calculation of the pre-filter 2 is included 
in the calculation of the transfer function identifying means 3. 

At first, transfer function to be identified in the first embodiment is set to be a quadratic model with making the wheel 
35 speed vibration Aco 1 detected by the pre-filter 2 at this time an output response, in the case of setting the road surface 
disturbance AT d to be the vibration input. Accordingly, the vibration model of the following formula is supposed. 



A*,--^ A7,+v, a 0 = l (15 ) 

i-0 

45 

In the above formula, v is an observed noise contained at a time of observing the wheel speed signal. When the formula 
(15) >s converted, the following formula can be obtained. 



50 z 2 

I^s^M^ = b 2 AT d + £a i s 2 "V (16) 

i-0 i-0 



55 

At first, after applying the pre-filter of the (14) formula to the (16) formula, the obtained formula is made discrete. At 
this time, Aco-i, AT d and v are expressed as the discrete data Aoo^kJ.ATcfk) and v(k) (in which k is a sampling number: 
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k = 1. 2, 3, ...) sampled at every sampling period T s . Further, the Laplace operator s can be made discrete by using a 
predetermined discrete method. In this embodiment, the discrete is performed by the following bilinear transformation 
as an example. In this case, d is one sample delay operator. 

2 l-d 

s = 



10 



15 



20 



25 



30 



35 



Further, since a dimension m of the pre-f ilter is desired to be two or more, in accordance with this embodiment, it 
is set that m = 2 by taking a calculating time into consideration, so that the following formula can be obtained. 



£a,£,tlc) = i> a &(*> + 5>i£i(*) <* 7 > 



2 

I 

i-0 i-0 

wherein 



= (S.) (1 + d) ' (1 -dr J r 0 (d)4mCk) (18) 

- ^ j (l + d)*F 0 (d)Ar d (k) (19) 

(l+d) i (l-d) J_1 F 0 (dV(*) (20) 

F 0 (d) = - ; — (21) 



Still further, in order to identify the transfer function from each of the data of the wheel speed vibration Aca-j on the 
basis of the least square method, the (1 7) formula is converted into the following formula in such a manner as to be a 
40 linear function type with respect to the parameter to be identified. In this case, " r is set to be a transposition of a matrix 

€yo{k) » C{k)6 ♦ r{k) (22) 
45 wherein 

8 - [an a 2 ] r 

SO * 

r(k) ^^W + Xa^.W (23) 

i-0 



55 

In the above formula, 0 is a parameter of a transfer function to be identified. 
Next, an operation of this embodiment will be described below. 

In the transfer function identifying means 3, the unknown parameter e is estimated by applying the least square 
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10 



35 



45 



50 



55 



method to each of the data obtained by successively applying the discrete data of the detected wheel speed vibration 
Ag>i to the formula (22), thereby identifying the transfer function. 

Concretely speaking, the detected wheel speed vibration A<d 1 is converted into the discrete data A<o (k) (k = 1 , 2, 3, 
...), the data is sampled at N points and the calculating formula of the least square method shown by the following for- 
mula is used, so that the parameter e of the transfer function is estimated. 

ks) - [t ff- k ak)C (*)J l |i p"- k C(*)£ 0 oc)J ( 24 ) 



15 Here, an amount to which a symbol A " is added is defined to be an estimated value thereof. 

Further, the least square method mentioned above may be calculated as a successive type least square method 
of calculating the parameter 0 by the following recurrence formula. 

20 * * 

6(N) - 6(N - 1) + h(N)[$M - CiNMN ~ l)] (25) 

25 P(N) - (J - h(N)C(N) \P(N-1) (21) 

P 



30 

In which p is a so-called forgetting factor, and is normally set to be a value of 0.95 to 0.99. At this time, an initial 
value may be set to be the following formula. 



- 0, P(-l)=ai. a: sufficiently large positive number 



40 Further, as a method of reducing an estimating error in the least square method mentioned above, various kinds of 
modified least square method may be employed. In this embodiment, an example of using the method of instrument 
variable corresponding to the least square method employing instrument variable will be described. In accordance with 
the method, with setting m(k) to be instrument variable at the stage of obtaining the relation of the (22) formula, the 
parameter of the transfer function is estimated by using the following formula. 



k*) - [l^'*in(Jc)rCk)J |t ^"*inCk)$ 0 Ck)J (28) 



Still further, the successive calculation is shown as follows. 
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6(N) = 6(N - 1) + h(N)[$ 0 (lt) - C(N)e(N - 1)] (29) 

h(lf)= P{N~1MN) (30) 

P+C(N)P(N -l)m(W) 

P(N) - A {J ~ h(N)f(N) \P(N - 1) (31) 
P 

JO 

A principle of the method of instrument variable is described as follows. Since the (22) formula is given to the (28) 
formula, the following formula is obtained. 

is 

W = 0 + |tp"~Vic)rOoJ |l^"'«C*)rtk)J (32) 



20 



the estimating value of 8 coincides with a true value of 0 when the observed value of concomitant variable is selected 
such that a second term of a right side in the (32) formula becomes zero. Then, in accordance with this embodiment, 
25 an observed value of concomitant variable delaying C(k) = R y1 (k) - 4 y2 (k)] T at a degree of having no relation with 
respect to a formula error r(k), is employed as the observed value of concomitant variable. 
Accordingly, it is set to be 



30 



35 



40 



m(k) = [^ y1 (k-L)-4y2(k-L)] T (33). 



In which L is a number of delay. 

As mentioned above, after identifying the transfer function, in the u gradient calculating means 4, the physical quan- 
tity in connection with the gradient D 0 of the road surface p. is calculated as follows. 



D * (34) 



a, J x +J 2 



When the physical quantity in connection with the gradient D 0 of the road surface n can be calculated by the (34) for- 
mula as mentioned above, for example, the frictional characteristic between the tire and the road surface can be easily 
diagnosed to be in a saturated state when the physical quantity is small. 

45 Next, a second aspect will be described below with reference to Fig. 9. As shown in Fig. 9, in accordance with the 

second aspect, adapting means 5 for changing the characteristic of the pre-filter 6 in accordance with the parameter 
identified by the transfer function identifying means 3 is further provided. 

In the second aspect, it is desirable that the transfer function of the pre-filter 6 has the same characteristic as that 
of a polynomial expression in a denominator of the transfer function to be identified. However, since the transfer function 

so to be identified is unknown, the transfer function of the pre-filter 6 is constituted by the following formula by using the 
parameter of the transfer function identified by the transfer function identifying means 3. 



55 
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F(s) = — 



i-0 



(35) 



10 

In the transfer function identifying means 3 at that time, the following calculation is performed in the same manner 
as the successive type least square method in accordance with the first aspect. 



15 



20 



30 



6(N) = d(N " 1) + h(N)[$ 0 (N) -CiNWN-l)] (36) 

W - (37) 

P^^ll-hWrWlPtW-l) (38) 

p 



wherein 

25 C(ie) - [-f y i(JO -f y2 (*)] r (39) 

^(Jk) = (1 + d) J (1 - d) 2 " 1 F 0 (d)^a) l (Jc) (40) 

P.(d)-- ^ 77^ (41) 



Then, in the adapting means 5, the coefficient of the transfer function of the (35) formula of the pre-filter 6 is con- 
40 structed from each of the components sl^ and a 2 of the parameter e of the estimated transfer function. Then, the transfer 
function identifying means 3 again estimates the parameter of the transfer function of the wheel resonance system from 
the wheel speed vibration Ao>i detected by the pre-filter 2 to which the coefficient is adapted. 

In accordance with the second aspect, since the characteristic of the pre-filter is properly changed in correspond- 
ence to the estimated parameter, there is an effect that a better estimated value can be obtained in comparison with the 
45 first aspect in which the parameter of the pre-filter is fixed. 

An estimating result in the case of performing the second aspect is shown in Fig. 10. In the same drawing, a hori- 
zontal axis indicates a true value of the physical quantity D 0 in connection with the easiness of the slip between the tire 
and the road surface, and a vertical axis indicates an estimated value of D 0 in accordance with the second aspect. As 
shown in the drawing, it is understood that the physical quantity D 0 in connection with the easiness of the slip between 
so the tire and the road surface is well estimated by the present aspect. 

In the case of applying the present invention to the antilock brake control apparatus, the Braking force is controlled 
such that the estimated physical quantity D 0 coincides with the reference value. Since the physical quantity D 0 is calcu- 
lated without premising that the resonance frequency is known, the antilock brake controlling operation can be accu- 
rately performed without relation whether or not the resonance frequency is changed by the replacement of the tire, the 
55 reduction of the tire air pressure or the like. Further, the same effect can be obtained also in the application to the trac- 
tion control. 
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(Second Embodiment) 



As is also understood from the (6) and (12) formulas, a constant K of the torsion spring of the tire is included in the 
vibration model in accordance with the present invention as the unknown component of the wheel state. With paying an 
5 attention to this point, in accordance with the second embodiment, in the same manner as the first embodiment it is 
structured such as to estimate the physical quantity in connection with the torsion spring constant K of the tire from the 
transfer function of the vibration model in the case that only the road surface disturbance AT d is input to the wheel res- 
onance system as the vibration input. 

The second embodiment will be described below with reference to a block diagram shown in Fig. 11 . In this case, 
io the same reference numerals are attached to the same components as those of the first embodiment, and a detailed 
explanation thereof will be omitted. 

As shown in Fig. 11, the second embodiment is structured such as to connect the spring constant calculating 
means 7 for calculating the physical quantity in connection with the torsion spring constant of the tire on the basis of the 
parameter identified by the transfer function identifying means 3 to the output end of the transfer function identifying 
is means 3. 

Next, a flow of the calculation by the spring constant calculating means 7 will be described below. 
When the parameter of the transfer function is identified by the transfer function identifying means 3 shown in Fig. 
1 1 , the identified transfer function can be expressed by the following formula. 



20 



G(s) = 



(42) 



25 Accordingly, from the corresponding relation to the vibration model expressed by the (1 2) formula, the following formula 
can be obtained. 



30 



a 2 =K/J l (43) 



Then, in the spring constant calculating means 7, in accordance with the (43) formula, the torsion spring constant K of 
the tire is calculated from the estimated value ( A is added) of the parameter a 2 identified by the transfer function identi- 
35 tying means 3 and the predetermined moment of inertia Jj in the rim end. 

Further, the following formula expresses the resonance frequency of the wheel resonance system. 

^= (44) 

40 

Then, in the spring constant calculating means 7, the resonance frequency of the wheel resonance system can also be 
calculated by the (44) formula as one of the physical quantities in connection with the torsion spring constant K of the 
tire. 

45 Further, since the torsion spring constant K of the tire and the resonance frequency obtained by this embodiment 
are dependent on the tire air pressure, they can be used as the information for a tire air pressure judging apparatus, so 
that a safe running can be performed. 

In accordance with this tire air pressure judging apparatus, for example, in the case that the estimated torsion 
spring constant K of the tire and the predetermined reference value are compared with each other and the constant is 

so equal to or more than the reference value, it is judged that "tire air pressure is normal", and in the case that the constant 
is less than the reference value, it is judged that "tire air pressure is abnormal", thereby warming the driver. Further, it 
may be structured such that a second reference value corresponding to the tire air pressure suitable for running on a 
high way is prepared, and in the case that the spring constant K is less than the second reference value, a warning of 
supplementing the tire air pressure is indicated to the driver. In this case, by preparing the reference values for various 

55 kinds of tires, the reference for judging the tire air pressure can be kept the same even at a time of replacing the tire. 

Further, there is an advantage that the vibration frequency can be provided to the antilock brake control apparatus 
for controlling the braking force on the basis of the vibration of the Braking force. For example, even when the resonance 
frequency changes due to the change of the tire air pressure and the replacement of the tire, the resonance gain cal- 
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culated in the same frictional state indicates substantiatty the same value by finely oscillating the Braking force by the 
changed resonance frequency, so that the detecting accuracy of the peak n can be maintained. Still further, in the anti- 
lock brake control apparatus for detecting a state immediately before the peak n by detecting the resonance frequency 
from the peak value of the frequency spectrum in the wheel speed signal and comparing the detected resonance fre- 

s quency with the resonance frequency when the tire grips, a state immediately before the peak \i can be always detected 
in an accurate manner without relation to the change of the tire air pressure or the like by setting the resonance fre- 
quency estimated by this embodiment to be a reference for comparison. 

In this case, in accordance with the second embodiment, the parameter of the pre-f ilter 2 is fixed. However, it may 
be structured such as to properly change the coefficient of the pre-f ilter 2 in accordance with the parameter of the esti- 

10 mated transfer function as in the second aspect of the first embodiment. 

An estimating result of the resonance frequency at a time of performing the second embodiment is shown in Fig. 
12. In the drawing, a horizontal axis indicates a true value of the resonance frequency, and a vertical axis indicates the 
resonance frequency calculated by using the (44) formula in accordance with this embodiment. As shown in the draw- 
ing, it is understood that the resonance frequency can be significantly well estimated by this embodiment. 

15 

(Third Embodiment) 

A wheel condition estimating apparatus in accordance with a third embodiment is structured such as to identify the 
transfer function of the wheel resonance system in the case that the vibration torque AT-j is input to the wheel reso- 
20 nance system as the vibration input, and has a first aspect for calculating the n gradient and a second aspect for calcu- 
lating the torsion spring constant of the tire. 

A structure of a first aspect of the wheel condition estimating apparatus in accordance with the third embodiment 
is shown in Fig. 13. In this case, the same reference numerals are attached to the same components as those of the 
first embodiment shown in Fig. 8, and the detailed explanation will be omitted. 
25 As shown in Fig. 1 3, a vibrating means 8 for adding the vibration torque AT 1 as the vibration input around the aver- 
age brake and drive torque T 1 is further provided in the wheel condition estimating apparatus. The vibrating means 8 
vibrates the wheel around the center axis thereof by vibrating the brake pressure and the engine output. For example, 
in the case of vibrating the Braking torque, the Braking torque can be vibrated by adding a pressure increase and 
decrease command having a predetermined frequency to an average brake pressure command in a control command 
30 of a pressure increase and decrease valve in a hydraulic actuator for a brake pressure. 

Next, a calculation formula on which the transfer function identifying means 3 in accordance with this embodiment 
depends will be introduced. In this case, also in accordance with this embodiment, the calculation of the pre- filter 2 is 
performed in such a manner as to be included in the calculation of the transfer function identifying means 3. 

At first, a quadratic model in which the transfer function to be identified in the third embodiment is made such that 
35 the vibration torque AT-, is a vibration input and the wheel speed vibration Ao>i detected by the pre-f ilter 2 at that time 
is an output response is prepared. Accordingly, the vibration model expressed by the following formula is supposed. 



2 



40 




(45) 



i-0 



45 

In which v is an observation noise included at a time of observing the wheel speed signal. The following formula is 
obtained by modifying the (45) formula. 

50 

Za^'yit) = £ b.s 2 " 1 ^ + £a iS s -V(t) (46) 

i-0 i-0 i-0 

55 



At first, after applying the pre-f ilter of the (14) formula to the (46) formula, the obtained formula is made discrete. At 
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this time, Ao^ , AT d and v are expressed as the discrete data Aa> 1 (k),AT d (k) and v(k) (in which k is a sampling number: 
k = 1 , 2. 3, ...) sampled at every sampling period T s . Further, the Laplace operator s can be made discrete by the bilinear 
transformation mentioned above using a delayed calculating component d. 

Further, since a dimension m of the pre-filter is desired to be two or more, in accordance with this embodiment, it 
5 is set that m = 2 by taking a calculating time into consideration, so that the following formula can be obtained. 



10 i-0 i-0 i-0 



(48) 



(49) 



(50) 
(51) 



15 

wherein 



20 
25 

F 0 (d) = 

30 



= [^"] (l + d?(X -d) r, F,(d)AK*) 

= {j^j (l + d) i (i -d) i_i F 0 (dM*) 
l 

i-o \ 2 / 



Still further, in order to identify the transfer function from each of the discrete data of the wheel speed vibration A© 1 
35 on the basis ol the least square method, the (47) formula is converted into the following formula in such a manner as to 
be a linear function type with respect to the parameter to be identified. In this case, " r is set to be a transposition of a 
matrix. 



40 

Zyo(k) - C(k)d + r(k) (52) 
wherein 

0 - [aia 2 ] T 

rOc) - ib&Qc) + i>x£i(*> < 53 > 

50 i-0 i-0 



In the above formula. 6 is a parameter of a transfer function to be identified. 
55 Next, an operation of this embodiment will be described below. 

In the transfer function identifying means 3, the unknown parameter 0 is estimated by applying the least square 
method to each of the data obtained by successively applying the discrete data of the detected wheel speed vibration 
AO! to the formula (52), thereby identifying the transfer function. 
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Concretely speaking, the detected wheel speed vibration Ao^ is converted into the discrete data Ao (k) (k = 1 , 2, 3, 
...). and the data is sampled at a number of N points. Then, by using the above formulas C(k) and £ y o(k), the estimated 
values < A is added) of the parameters and a 2 of the transfer function are calculated in accordance with the same cal- 
culation as that of the (24) formula and thereafter in the first embodiment. In this case, in this embodiment, the succes- 
s sive type least square method and the method of observed value of concomitant variable may be employed, and as in 
the second aspect of the first embodiment, the coefficient of the pre-f ilter 2 may be properly changed in correspondence 
to the parameter of the identified transfer function. 

In this case, an example of the signal waveform of the vibration torque AT 1 given by the vibrating means 8 will be 
shown in Figs. 14A, 14B and 14C. 
10 The signal in Fig. 1 4A is a pseudo random signal, which, for example, can be generated on the basis of an M group 
signal as is well known as one of the pseudo random signals. In this case, since the vibration input has a lot of frequency 
components, there is an advantage that an accuracy of estimating the transfer function is improved. 

Further, the signal in Fig. 14B is a step signal which suddenly stands up at a certain time from a state in which the 
vibration torque is 0 and keeps a constant vibration torque thereafter. In the case of using the step signal, since it is not 
15 necessary to frequently change the torque as in the pseudo random signal, the vibration and the noise are reduced, so 
that there is an advantage that the load of the actuator is reduced. 

Still further, the signal in Fig. 14C is an impulse signal which suddenly stands up at a certain time from a state in 
which the vibration torque is 0 and again returns to a state in which the vibration torque is 0 after a fixed time has 
passed. In the case of using the impulse signal, since there is no offset component in the torque and the signal can be 
20 immediately given in various running state such as the braking, the driving, the steady running and the like, the transfer 
function can be estimated in an optional running state. 

Next, as a second aspect, a wheel condition estimating apparatus further comprising the physical quantity in con- 
nection with the torsion spring constant of the tire in the vibration model in accordance with the first aspect mentioned 
above will be described below. 

25 As is understood from the formulas (5) and (1 1), in the vibration model of the present invention, the torsion spring 
constant K of the tire is included as an unknown component of the wheel condition. With paying attention to this point, 
in accordance with the second aspect, it is structured such as to estimate the physical quantity in connection with the 
torsion spring constant K of the tire from the transfer function of the vibration model in the case that the vibration torque 
AT! is input to the wheel resonance system as the vibration input. 

30 The second aspect will be described below with reference to a block diagram shown in Fig. 15. 

As shown in Fig. 15, the second aspect is structured such that in the same manner as that of the second embodi- 
ment, spring constant calculating means 7 for calculating the physical quantity in connection with the torsion spring con- 
stant of the tire on the basis of the parameter identified by the transfer function identifying means 3 is connected to an 
output end of the transfer function identifying means 3. 

35 Next, a flow of the calculation by the spring constant calculating means 7 will be described below. 

When the parameter of the transfer function is identified by the transfer function identifying means 3 shown in Fig. 
1 5, the identified transfer function can be expressed by the following formula. 

4 °- G(s) s a » (54) 

s + a x s + a 2 



Accordingly, the following formula can be obtained from the corresponding relation to the vibration model expressed by 
45 the (11) formula. 

a 2 =KlJ x (55) 

50 

Then, in the spring constant calculating means 7, in accordance with the (55) formula, the torsion spring constant K of 
the tire is calculated from the estimated value ( A is added) of the parameter a 2 identified by the transfer function identi- 
fying means 3 and the predetermined moment of inertia Ji in the rim end. 

Further, the following formula expresses the resonance frequency of the wheel resonance system. 
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yfa^= Jk/J x (56) 

5 

Then, in the spring constant calculating means 7, the resonance frequency of the wheel resonance system can also be 
calculated by the (44) formula as one of the physical quantities in connection with the torsion spring constant K of the 
tire. 

Further, since the torsion spring constant K of the tire and the resonance frequency obtained by this aspect are 
w dependent on the tire air pressure, they can be used as the information for a tire air pressure judging apparatus. Fur- 
ther, there is an advantage that the vibration frequency can be provided with respect to the antilock brake control appa- 
ratus for controlling the braking force on the basis of the vibration of the braking force. 

In this case, in accordance with this aspect, the signal waveform of the vibration input can be made any one of the 
waveforms shown in Figs. 14A, 14B and 14C as in the first aspect. 
is As mentioned above, in the third embodiment, since the output response equal to or more than a fixed value can 
always be obtained by applying the vibration torque AT 1 to the wheel resonance system even in the case that the ampli- 
tude of the wheel speed vibration Aa>i is made small due to the little unevenness of the road surface, the estimation of 
the n Gradient and the torsion spring constant of the tire can be accurately performed without relation to the road sur- 
face condition. 

20 

(Fourth Embodiment) 

A wheel condition estimating apparatus in accordance with a fourth embodiment is structured such as to identify 
the transfer function of the wheel resonance system from the detected vibration input and output response in the case 

25 that the vibration torque AT 1 is input to the wheel resonance system as the vibration input, and has a first aspect for 
calculating the \x Gradient and a second aspect for calculating the torsion spring constant of the tire. 

A structure of a first aspect of the wheel condition estimating apparatus in accordance with the first aspect of the: 
fourth embodiment is shown in Fig. 16. In this case, the same reference numerals are attached to the same compo- 
nents as those of the third embodiment shown in Fig. 13, and the detailed explanation will be omitted. 

30 As shown in Fig. 16, this wheel condition estimating apparatus is further provided with a vibration input detecting 
means 9 for detecting an actual value of the vibration torque AT^ as the vibration input given to the wheel by the vibrat- 
ing means 8. Then, the pre-filter 2 gives a predetermined filter treatment to the detected vibration input together with 
the wheel speed , and the transfer function identifying means 3 estimates the parameter of the transfer function from 
the filter-treated vibration input component and the wheel speed vibration Ao^ 

35 A method of detecting the vibration input by the vibration input detecting means 9 is structured, for example, in the 
case of vibrating the Braking torque, such as to detect the brake pressure (the wheel cylinder pressure) corresponding 
to the Braking torque by means of the pressure sensor or the like, to convert it into the Braking torque including the 
vibration component by multiplying a predetermined constant, and to subtract the value of the average Braking torque 
from this Braking torque, so that only the vibration torque component is extracted as the vibration input. 

40 Next, a calculation formula on which the transfer function identifying means 3 in accordance with this embodiment 
depends will be introduced. In this case, also in accordance with this embodiment, the calculation of the pre-filter 2 is 
performed in such a manner as to be included in the calculation of the transfer function identifying means 3. 

At first, a quadratic model in which the transfer function to be identified in the fourth embodiment is made such that 
the vibration torque AT 1 is a vibration input and the wheel speed vibration Ao^ detected by the pre-filter 2 at that time 

45 is an output response is prepared. Accordingly, the vibration model expressed by the following formula is supposed. 



Ao\ = *f 4r l +v,a 1) =l (57) 

i-0 



55 

In the above formula, AT-, is an actual vibration torque detected by the vibration input detecting means 9 and v is an 
observation noise included at a time of observing the wheel speed signal. The following formula is obtained by modify- 
ing the (57) formula. 
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i-0 i-0 



At first, by applying the pre-filter of the (14) formula to the (58) formula, the obtained formula is made discrete. At 
10 this time, Ao^, AT d and v are expressed as the discrete data Aco^kJ.ATdtk) and v(k) (in which k is a sampling number: 
k m i , 2, 3, ...) sampled at every sampling period T s . Further, the Laplace operator s can be made discrete by the bilinear 
transformation mentioned above using a delayed calculating component d. 

Further, since a dimension m of the pre-filter is desired to be two or more, in accordance with this embodiment, it 
is set that m = 2 by taking a calculating time into consideration, so that the following formula can be obtained. 

75 



£«t£i(k) " 2><U*) + (59) 



20 



25 



30 



35 



2 

I 

i-0 1-0 i-0 

wherein 



COc) - fel (1 + d) 1 (l - d) J " i F 0 (d)^T 1 (k) (61) 
=^J(l + d) i (l-d) l - i P 0 (d) V (k) (62) 
- - n , 1 (63) 

£ Cl p- (l+dWi-d)*" 1 

i-o \ 2 / 



Still further, in order to identify the transfer function from each of the discrete data of the wheel speed vibration Aa>i 
40 on the basis of the least square method, the (47) formula is converted into the following formula in such a manner as to 
be a linear function type with respect to the parameter to be identified. In this case, " T ' is set to be a transposition of a 
matrix. 



? y o - f(k)6 + r{k) (64) 
wherein 

so C(*> - l-fwU> '€y2lk){ u o{k)€ al (k)€ u2 ik)} T 

9 = [ai a2 bo bi b2] T 

r(k) = L^CiOc) (6 5) 

55 
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In the above formula, 0 is a parameter of a transfer function to be identified. However, in accordance with this embodi- 
ment, since the AT 1 is detected, it is understood that the coefficients b 0 , b-i and b 2 can be also estimated as is different 
from the embodiment mentioned above. 

Next, an operation of a first aspect in accordance with the fourth embodiment will be described below. 

5 In the transfer function identifying means 3, the unknown parameter 6 is estimated by applying the least square 

method to each of the data obtained by successively applying the discrete data of the detected wheel speed vibration 
Acot to the formula (64), thereby identifying the transfer function. 

Concretely speaking, the detected wheel speed vibration Ag>^ is converted into the discrete data Aa>(k) (k = 1 , 2, 3, 
...), and the data is sampled at a number of N points. Then, by using the above formulas £(k) and 4 y0 (k). the estimated 

w values ( A is added) of the parameters and a 2 of the transfer function are calculated in accordance with the same cal- 
culation as that of the (24) formula and thereafter in the first embodiment. In this case, in this embodiment, the succes- 
sive type least square method and the method of observed value of concomitant variable may be emptoyed. In the case 
of using the method of observed value of concomitant variable, when the vibration input can be detected as in this 
embodiment, the estimated value ( A is added) of the true output x(k) is successively calculated by using the parameter 

is of the identified transfer function in accordance with the following formula 



20 x(k) ■ '~* ) 2 £ ^ (66) 

i-o \2J 

25 

Accordingly, the observed value of concomitant variable can be constructed in the following formula by using this esti- 
mated value. 



30 y 

• u*) = ci+dyft-dr^d^ck) (67) 

*(*) » [-&!<*) -? x *(k) £ u0 ik) t vl (k) £ u2 (k)] T (68) 

35 



40 Further, a formula structured such as to delay the following formula in such a manner as to have no relative relation 
with the formula error r(k) is used. 

45 Accordingly, the following formula is given, in which L is a delayed time. 

m(k) = [- !, y1 (k - L) - t y2 (k - L) ^ u0 (k - L)^ u1 (k - LK u2 {k - L)] 7 (70) 

Further, in accordance with this aspect, as in the second aspect of the first embodiment, the coefficient of the pre- 
50 filter 2 can be properly changed in correspondence to the parameter of the identified transfer function. 

Still further, the waveform of the vibration input may be structured in such a manner as shown in Figs. 14A, 14B 
and 14C. In the case that the resonance frequency at a time when the tire is gripped on the road surface in this aspect 
(hereinafter, refer to "tire resonance frequency") is apparently known, it may be vibrated by the tire resonance fre- 
quency. In accordance with the experimental research by the inventors, it is known that the vibration model expressed 
55 by the formula (1 1 ) becomes the following formula. 

J, - J 2 (71) 
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Accordingly, when the transmitting characteristic from the vibration input to the output response at a time of being 
vibrated by the tire resonance frequency is calculated on the basis of the can be obtained. 



Aw = -J^l-A^ (72) 
2KJ X 

10 

Then, the transfer function is identified as the following formula by using the least square method in the transfer function 
identifying means 3 by setting the transfer function to be identified as Aco 1 = G d AT 1 (73) . 



15 



20 



25 



30 



6(n) = [i^CCWrcwj |£^"*C0fc)$ o (k)J (74) 

wherein 

e = G, (75) 

£ 0 (Jc) - Awjk) (7 6) 

&Cfc)-4qCk) (77) 



Further, the least square method mentioned above can be calculated as the successive type least square method 
35 as in the following formula. 

- 0(N - 1) + h(W)[< 0 (w) ~ C (N)d(N " 1)] .(78) 



40 



45 



P(N) = [I - h(N)C(N) )P{N - 1) ( 80 ) 

P 



so In the above formula, p is a so-called forgetting factor, 

In the above formula, p is a so-called forgetting factor, and is normally set to be a value of 0.95 to 0.99. At this time, 
an initial value may be set to be the following formula. 

55 ~ 0, P(-l)=cci, a: sufficiently large positive number 

Further, a further accurate estimated value can be obtained by applying a modified least square method such as the 
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method of observed value of concomitant variable or the like to the least square method mentioned above. 

Then, the estimated value of the parameter Gd of the identified transfer function is expressed by the following for- 
mula in accordance with the corresponding relation between the formulas (72) and (73). 




(81) 



10 Accordingly, the \i gradient calculating means 4 calculates the physical quantity in connection with the easiness of the 
slip D 0 of the road surface in accordance with the (81 ) formula by using the estimated value of the Gd. 

In this case, the frequency of the vibration input may be the resonance frequency calculated by the torsion spring 
constant of the tire which is calculated in accordance with the second embodiment of the present invention and the sec- 
ond aspect of the third embodiment. In the case of vibrating by only the tire resonance frequency as mentioned above, 
is since the power of the vibration input can be concentrated to the frequency component, an SN ratio of the output from 
the vibration input detecting means 9 and the output response detecting means 1 is improved, so that there is an advan- 
tage that the calculating accuracy is improved. 

Next, as a second aspect, a wheel condition estimating apparatus further comprising a physical quantity in connec- 
tion with the torsion spring constant of the tire in the vibration model in the first aspect mentioned above, will be 
20 described below with reference to a block diagram shown in Fig. 17. 

As shown in Fig. 17, the second aspect is structured such that in place of the \x gradient calculating means 4 in 
accordance with the first aspect, the spring constant calculating means 7 for calculating the physical quantity in connec- 
tion with the torsion spring constant of the tire on the basis of the parameter identified by the transfer function identifying 
means 3 is connected to the output end of the transfer function identifying means 3. 
25 Next, a flow of the calculation by the spring constant calculating means 7 will be described below. 

When the parameter of the transfer function is identified by the transfer function identifying means 3 shown in Fig. 
1 7, the identified transfer function can be expressed by the following formula. 



Au x = *f +v,ao = 1 (82) 

35 i-0 



Accordingly, the following formula can be obtained from the corresponding relation to the vibration model expressed by 
the (11) formula. 

40 

^ = K (83) 

45 Then, in the spring constant calculating means 7, in accordance with the (83) formula, the torsion spring constant K of 
the tire is calculated from the estimated value ( A is added) of the parameter a 2 identified by the transfer function identi- 
fying means 3 and the predetermined moment of inertia in the rim end. 

Further, the following formula expresses the resonance frequency of the wheel resonance system. 

so 

Ja 2 = ylK/J t (84) 



55 

Then, in the spring constant calculating means 7, the resonance frequency of the wheel resonance system can also be 
calculated by the (44) formula as one of the physical quantities in connection with the torsion spring constant K of the 
tire. 
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Further, since the torsion spring constant K of the tire and the resonance frequency obtained by this aspect are 
dependent on the tire air pressure, they can be used as the information for a tire air pressure judging apparatus. Fur- 
ther, there is an advantage that the vibration frequency can be provided with respect to the antilock brake control appa- 
ratus for controlling the breaking pressure on the basis of the vibration of the Braking force. 

In this case, also in accordance with this aspect, the signal waveform of the vibration input can be made any one 
of the waveforms shown in Figs. 14A, 14B and 14C as in the first aspect except the case of vibrating by the tire reso- 
nance frequency. 

As mentioned above, in the fourth embodiment, since the vibration input actually detected is used, more transfer 
function parameters can be estimated. 

(Fifth Embodiment) 

A wheel condition estimating apparatus in accordance with a fifth embodiment is structured such as to select an 
output response including a large amount of a periodical signal among the output responses, identify a transfer function 
in a wheel resonance system on the basis of the selected output response, and calculate a \x gradient. 

A structure of the wheel condition estimating apparatus in accordance with the fifth embodiment will be described 
below with reference to a block diagram shown in Fig. 18. As shown in Fig. 18, the structure is the same with the struc- 
ture of the first embodiment except that selecting means 1 0 is provided between the pre-f ilter 2 and the transfer function 
identifying means 3, the selecting means 10 selects a periodical signal among the wheel speed vibration Ao-j output 
from the pre-f ilter 2 and outputs the selected output response to the transfer function identifying means 3. 

In the case that the wheel speed vibration Acd 1 is a periodical signal, when the wheel speed vibrations Aco 1 input for 
a predetermined period are plotted on the frequency-amplitude (power) graph, it is expected that a graph shown by a 
reference numeral 100 in Fig. 19 can be obtained. 

However, When the transfer function in the wheel resonance system is identified by using an output 102 having an 
unexpectedly large amplitude (for example, a returning noise by an aliase and the like) and an output 104 in which a 
periodical change does not significantly appear (an output having a small S/N ratio) ,and the p gradient is calculated, 
there is a great dispersion and an estimating accuracy is deteriorated. Accordingly, it is necessary to select the output 
response corresponding to the periodical signal. 

In comparison with Fig. 20B which graphs the periodical output 100 in a time sequence manner, the output 102 
having an unexpectedly large amplitude has a bad periodical characteristic when graphed in a time sequence manner 
as shown in Fig. 20A. Further, the output 104 in which the periodical change does not significantly appear has a very 
bad periodical characteristic when graphed in a time sequence manner as shown in Fig. 20C. 

On the contrary, the periodical output 100 shown by a solid line in Fig. 20B and the value (shown by a dot line in 
Fig. 20B) obtained by converting the periodical output 100 such that the periodical portion of the periodical output 100 
changes in correspondence to the change of the periodical output 100 have a certain relation. 

Accordingly, as explained in the embodiment mentioned above, when the periodical output 100 is set to be £ y0 , and 
the value obtained by converting the periodical output 100 in the above manner, for example, performing an even 
number of (for example, two times) differentiation (or any of an even number of difference, an even number of integration 
and an even number of summation) is set to be ^ a point defined by ^ (= x) and ^ (= y) is positioned within a cer- 
tain range 1 1 0 as shown in Fig. 21 . When the output is completely periodical, the point is positioned on a predetermined 
straight line 1 12 and the point is positioned in a wider range 1 14 as the periodical characteristic becomes worse. 

Then, in accordance with this embodiment, the apparatus is structured such as to calculate the relational coeffi- 
cient y (the following formula (85)) between the output response £yo corresponding to the wheel speed vibration Aoo-i 
and the value ^ obtained by converting the output response ^ in a manner mentioned above, and select the output 
response corresponding to the periodical signal on the basis of the relational coefficient a. 



a - 




(85) 



In the above formula, p c is a forgetting factor and N is a number of the output responses. 



24 



.0891904A2J_> 



EP 0 891 904 A2 



Accordingly, the selecting means 1 0 performs a selecting operation routine shown in Fig. 22 at every time when the 
wheel speed vibration A<o 1 is input from the pre-lilter 2. 

In a step 22, a pre-filter operation (a band area passing operation) is performed. In this case, a band area of the 
pre-filter operation is wider than a band area of the pre-filter 2. 

5 In a step 24, the relational coefficient a (the formula 185)) between the output response ^ corresponding to the 

wheel speed vibration Aa> 1 after the pre-filter operation and the value %p obtained by converting the output response 
£yo in a manner mentioned above is calculated. 

In a step 26, whether or not the wheel speed vibration Aco-j input at this time is a periodical signal is determined by 
determining whether or not the calculated relative coefficient a is within a predetermined range (a t < a < a 2 ). 

w In the case that the calculated relative coefficient a is not within the predetermined range, the wheel speed vibration 
A©! input at this time is not a periodical signal, so that this routine is finished. Accordingly, since the wheel speed vibra- 
tion A©-, input at this time is not input to the transfer function identifying means 3, it is prevented that the transfer function 
in the wheel resonance system is identified on the basis of the non-periodical signal so as to calculate the \i gradient. 
On the contrary, in the case that the calculated relative coefficient a is within the predetermined range, the wheel 

is speed vibration A©! input at this time is a periodical signal, so that in a step 28, the wheel speed vibration Aco 1 input at 
this time is input to the transfer function identifying means 3. Accordingly, the transfer function of the wheel resonance 
system is identified on the basis of the periodical signal and the *i gradient is calculated. Therefore, an accuracy of iden- 
tifying the transfer function and an accuracy of estimating the \x gradient can be improved. 

In this case, even when the output having an unexpectedly large amplitude (a noise) is converted in the above man- 

20 ner, the converted value includes a lot of components other than the periodical portions. Accordingly, since the value 
obtained by converting the output having an unexpectedly large amplitude in the manner mentioned above does not 
change in correspondence to the change of the output response, the relative coefficient a between the output having 
an unexpectedly large amplitude and the value obtained by converting the output response in the manner mentioned 
above is not within the predetermined range. For example, as shown in Figs. 24A and 24B, a group of signals A(refer 

25 to Fig. 23A) having a frequency including a noise K within the predetermined range containing a frequency and hav- 
ing a relative coefficient a outside a predetermined range is removed by the selecting means 10. 

Further, as shown in Fig. 25. in the value obtained by converting the output C in which a periodical change does 
not significantly appear in the manner mentioned above, there is a small number of components converted such as to 
change in correspondence to the change of the output response. Accordingly, the correlation coefficient a between the 

30 output C in which the periodical change does not significantly appear and the converted value obtained by converting 
the output response in the manner mentioned above does not fall in the predetermined range mentioned above. For 
example, as shown in Figs. 26A and 26B, a group of signals C having correlation coefficient a falling outside the prede- 
termined range (see Fig. 26A) is removed by the selecting means 10. The output C in which the periodical change does 
not significantly appear was obtained when the vehicle ran on a dirt road. 

35 On the contrary, as shown in Figs. 24C and 24D, a group of signals B (refer to Fig. 23A) having no noise component 
and having a relative coefficient a within a predetermined range is selected by the selecting means 1 0 so as to be input 
to the transfer function identifying means 3. As shown in Fig. 25, in the value obtained by converting the periodical out- 
put D in the manner mentioned above, there is a large number of components converted such as to change in corre- 
spondence to the change of the output response. Accordingly, the correlation coefficient a between the periodical 

40 output D and the value obtained by converting the output response in the manner mentioned above falls in the prede- 
termined range mentioned above. For example, as shown in Figs. 26A and 26B, a group of signals D having correlation 
coefficient a falling within the predetermined range (see Fig. 26A) is not removed by the selecting means 10. Accord- 
ingly, the group of the signals D are input to the transfer function identifying means 3. The periodical output D shown in 
Fig. 25 was obtained when the vehicle ran on an asphalt road. 

45 Accordingly, the unknown component mentioned above can be estimated on the basis of the output response from 
which the output response having an unexpectedly large magnitude and the output in which the periodical change does 
not significantly appear are removed. 

In accordance with the fifth embodiment mentioned above, the \i gradient is calculated. However, the present inven- 
tion is not limited to this, and a torsion spring constant of a tire may be calculated. Accordingly, this embodiment is not 

so applied only to the first embodiment (the first aspect) mentioned above, but can be applied to all the embodiments men- 
tioned above. 

In each of the embodiments in accordance with the present invention mentioned above, since the structure is made 
such as to identify the transfer function of the continuous time model and to directly estimate the parameter of the trans- 
fer function, the problems mentioned above which is generated together with identifying the discrete time model can be 
55 avoided. 

Further, since the wheel condition included in the vibration model established in the various kinds of running states 
is identified, the wheel condition in the various running conditions such as the braking, the driving, the steady running 
and the like can be estimated without limiting to the braking time. 
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The embodiments of the present invention are described above, however, the present invention is not limited to the 
embodiments mentioned above, and can be optionally and preferably changed within the scope of the present inven- 
tion. 

For example, in the first to fifth embodiments mentioned above, the explanation is given on the basis that the wheel 
5 resonance system is the quadratic vibration model expressed by the (1 1 ) and (12) formulas. However, it is apparent that 
the wheel condition estimation can be performed on the basis of the cubic vibration model expressed by the (5) and (6) 
formulas- 
Claims 

10 

1. A wheel condition estimating apparatus for estimating a wheel condition in a wheel resonance system including a 
frictional characteristic between a tire and a road surface, comprising: 

detecting means for detecting an output response with respect to a vibration input to the wheel resonance sys- 
15 tern; and 

estimating means for expressing a transfer characteristic of said wheel resonance system from the vibration 
input to the output response by a vibration model including, at least a physical quantity relating to ease of a slip- 
ping between the tire and the road surface as an unknown component of a wheel condition, and for estimating 
at least said unknown component which substantially satisfies the output response detected by said detecting 
20 means on the basis of the vibration model. 

2. A wheel condition estimating apparatus according to claim 1 , wherein said model is a transfer function of said 
wheel resonance system from the vibration input to the output response, and said estimating means estimates said 
unknown component by identifying said transfer function. 

25 

3. A wheel condition estimating apparatus according to claim 2, wherein said estimating means expresses said trans- 
fer function by a function of a parameter for expressing said physical quantity as said unknown component, esti- 
mates said parameter on the basis of said function and the output response detected by said detecting means, and 
estimates said unknown component on the basis of the estimated parameter. 

30 

4. A wheel condition estimating apparatus according to claim 1, further comprising modifying means for modifying 
said vibration model on the basis of the output response detected by said detecting means. 

5. A wheel condition estimating apparatus according to claim 1 , wherein said vibration model further includes, a phys- 
35 ical quantity relating to a constant of a torsion spring of the tire as an unknown component of the wheel condition. 

6. A wheel condition estimating apparatus according to claim 1, further comprising vibrating means for applying a 
vibration input to said wheel resonance system. 

40 7. A wheel condition estimating apparatus according to claim 6, further comprising vibration input detecting means for 
detecting the vibration input applied to said wheel resonance system by said vibrating means. 

8. A wheel condition estimating apparatus according to claim 1 , wherein said vibration input is at least one of a vibra- 
tion torque to said wheel resonance system and a road surface vibration torque from a road surface to said wheel 

45 resonance system. 

9. A wheel condition estimating apparatus according to claim 1, wherein said output response is a wheel speed. 

10. A wheel condition estimating apparatus according to claim 1, further comprising: 

50 

calculating means for calculating a determining value for determining whether or not the output response 
detected by said detecting means is a periodical output on the basis of said output response; and 
selecting means for selecting an output response to be input to said estimating means among the output 
responses detected by said detecting means on the basis of the determining value calculated by said calculat- 
55 ing means, wherein 

said estimating means uses the output response selected by said selecting means at a time of estimating said 
unknown component. 
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11. A wheel condition estimating apparatus according to claim 10, wherein said calculating means calculates a corre- 
lation coefficient between the output response detected by said detecting means and a converted value obtained 
by converting said output response such that the periodical portion of said output response changes in correspond- 
ence to the change of said output response, as said determining value. 

5 

12. A wheel condition estimating apparatus according to claim 1 1 , wherein said conversion is performed by any of an 
even number of differentiation, an even number of difference, an even number of integration and an even number 
of summation. 

10 1 3. A wheel condition estimating apparatus according to claim 1 2, wherein said calculating means calculates the rela- 
tive coefficient a from the following formula. 



15 
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20 



wherein £ y0 is the output response, £ y2 is the converted value, p c is a forgetting factor, N is a number of the output 
responses detected by said detecting means and a is the relational coefficient. 

25 14. A wheel condition estimating apparatus according to claim 10, wherein said calculating means calculates said 
determining value at every time when a predetermined number of output responses are detected by said detecting 
means. 

15. A wheel condition estimating apparatus according to claim 10, further comprising modifying means for modifying 
30 said vibration model on the basis of the output response detected by said detecting means. 

16. A wheel condition estimating apparatus according to claim 10, wherein said vibration model further includes, a 
physical quantity relating to a constant of a torsion spring of the tire as an unknown component of the wheel con- 
dition,. 

35 

17. A wheel condition estimating apparatus according to claim 10, further comprising vibrating means for applying a 
vibration input to said wheel resonance system. 

18. A wheel condition estimating apparatus according to claim 17, further comprising vibration input detecting means 
40 for detecting the vibration input applied to said wheel resonance system by said vibrating means. 

1 9. A wheel condition estimating apparatus according to claim 10, wherein said vibration input is at least one of a vibra- 
tion torque to said wheel resonance system and a road surface vibration torque from a road surface to said wheel 
resonance system. 



45 



50 



20. A wheel condition estimating apparatus according to claim 9, wherein said output response is a wheel speed. 
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